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Executive  Summary 


1 .  Objectives 

Techniques  for  controlling  the  behavior  of  turbulent  fluid  flows  are  beginning  to  be 
developed  as  a  result  of  fundamental  understanding  of  the  basic  mechanisms  active  in 
turbulent  flows.  For  example,  by  properly  combining  axial  and  helical  (orbital)  disturbances, 
one  can  cause  a  round  turbulent  jet  to  split  into  two  distinct  jets  ( bifurcating  jets )  or  to 
explode  into  a  shower  of  vortex  rings  forming  a  rapidly-mixing,  wide-angle  conical 
spreading  flow  ( blooming  jets).  These  phenomena  were  discovered  in  an  NSF-funded 
program  under  the  Principal  Investigator,  which  led  to  extensive  research  on  dual-mode 
forcing  under  AFOSR  support.  This  report  summarizes  the  research  accomplished  under 
AF-F49620-84-K-0005. 

The  general  objective  of  this  work  is  to  control  jet  flows  by  properly  organizing  the 
dominant  large-scale  vortex  structures.  The  work  has  focused  on  understanding  the 
influence  of  excitation  frequency  and  amplitude  on  the  large-scale  structures  of  these  flows 
at  various  Reynolds  numbers.  Specifically,  the  effects  of  dual-mode  forcing  on  the  structure, 
momentum,  and  mixing  characteristics  of  turbulent  jets  have  been  investigated  using  a 
combination  of  laboratory  experiments  and  numerical  simulations.  -  • 
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2.  Overview  of  results 

Dramatic  changes  in  jet  development  described  above  were  discovered  by  Lee  and 
Reynolds  (1985a;  see  Appendix  A)  in  a  study  of  mechanically-excited  water  jets  at  a 
Reynolds  number  of  4,000.  The  axial  excitation  was  provided  by  weak  pulsations  in  the  jet 
flow  and  the  orbital  excitation  by  mechanical  motion  of  the  nozzle.  The  basic  flow  regimes 
and  flow  structure  were  determined  using  flow  visualization  with  dye  and  laser-induced- 
fluorescence.  A  fiber-optic  laser-doppler-anemometry  system  was  used  to  measure  the 
velocity  field  and  momentum  mixing,  and  acid-suppressed  laser-induced  fluorescence  was 
used  to  study  molecular  mixing.  This  work  is  reported  in  detail  by  Lee  and  Reynolds 
(1985b).  ' 


To  provide  insights  into  the  governing  mechanisms  and  to  study  the  effect  of  various 
parameters,  a  vortex-filament  code  for  simulating  the  large-scale  characteristics  of  excited 
jets  was  developed  in  cooperation  with  NASA/Ames  (Leonard,  Couet,  and  Parekh  1985;  see 
Appendix  B).  This  work  provided  confirmation  of  the  mechanism  of  bifurcation  and 
blooming.  Concentric  vortex  rings,  the  dominant  flow  structures  in  the  near  field  of  a  normal 
jet,  leapfrog  over  one  another  axisymmetrically.  These  rings  are  made  eccentric  from  one 
another  by  orbital  excitation,  and  eccentric  rings  tend  to  tilt  one  another.  Under  the  right 
conditions,  this  interaction  forces  them  apart,  fit  dual-mode  excitation,  the  axial  excitation 
controls  the  phasing  of  the  vortex  rings,  and  the  orbital  excitation  controls  the  azimuthal 
phasing  of  the  rings.  If  two  rings  are  shed  per  orbit,  and  the  spacing  between  them  is  in  a 
certain  range,  the  jet  will  bifurcate.  If  the  frequency  ratio  is  a  fraction  in  a  certain  range,  the 
jet  will  bloom  into  a  shower  of  independent  vortex  rings. 

While  on  sabbatical  leave  at  Caltech  in  1984,  the  Principal  Investigator  constructed  a 
second  water  apparatus  in  which  bifurcation  and  blooming  were  demonstrated  with  a  smaller 
diameter  jet  that  allowed  visualization  up  to  100  diameters  downstream.  There  it  was  shown 
that  the  bifurcation  is  permanent,  and  hence  the  far  field  of  a  bifurcating  jet  is  dramatically 
different  than  that  for  a  normal  jet  with  the  same  momentum  flux. 
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The  experimental  results  in  the  water  facility  and  the  numerical  simulation  suggested 
that  similar  control  of  jet  flows  could  be  achieved  in  air  jets  at  higher  Reynolds  numbers. 
The  challenge  at  high  Reynolds  numbers  is  to  force  "collective  interactions"  of  the  closely- 
spaced  vortices  in  the  thin  shear  layer  occurring  at  high  Reynolds  number,  so  that  the  spacing 
of  the  resulting  large  vortices  is  placed  in  the  range  necessary  for  bifurcation  or  blooming.  A 
simple  air  flow  apparatus  was  constructed  in  which  it  was  found  that  weak  acoustic 
excitations  not  requiring  mechanical  nozzle  motion  would  produce  the  bifurcating/blooming 
phenomena  at  Reynolds  numbers  up  to  20,000.  (Parekh,  Reynolds,  and  Mungal  1987; 
Appendix  C).  Subsequently  another  air  apparatus  for  higher-speed  experimentation  was 
constructed,  and  it  was  found  that  the  phenomena  could  be  produced  by  acoustic  excitation  at 
a  Reynolds  number  of  100,000,  M  =  0.22,  where  the  vortex  rings  undergo  transition  to 
turbulence  before  the  bifurcation  or  blooming  (Parekh  and  Reynolds  1987;  Appendix  D). 
This  provided  the  first  clear  evidence  that  the  phenomena  could  be  effected  in  flows  of 
primary  technical  interest  to  the  AFOSR. 

3.  Summary  of  principal  conclusions 

1.  Properly  combined  axial  and  helical  excitations  can  produce  significant  changes 
in  the  evolution  of  round  turbulent  jets,  including  spreading  angles  as  large  as  80 
degrees. 

2.  Bifurcating  and  blooming  jets  are  dominated  by  vortex  rings  which  evolve  in 
well-defined  patterns  governed  by  the  ratio  of  axial-to-helical  excitation 
frequencies. 

3.  Bifurcating  jets  occur  when  the  axial/orbital  frequency  ratio  is  2  and  the  axial 
Strouhal  number  is  between  0.3  and  0.7.  Within  this  range,  the  spreading  angle 
increases  with  axial  Strouhal  number. 

4.  The  two  branches  of  the  bifurcating  jet  evolve  as  separate  jets,  apparently 
indefinitely  far  from  the  jet  exit. 

5.  Blooming  jets  occur  when  the  axial/orbital  frequency  ratio  is  non-integer 
between  1.7  and  2.8  and  the  axial  Strouhal  number  is  between  0.3  and  0.7. 
Within  this  range,  the  spreading  angle  increases  with  axial  Strouhal  number. 

6.  Molecular  mixing  is  enhanced  in  blooming  jets. 

7.  Increasing  either  the  axial  or  helical  excitation  amplitude  can  cause  the  spreading 
angle  of  the  bifurcating  jet  to  increase.  However,  increasing  amplitudes  too  high 
results  in  saturation  that  inhibits  further  increases  in  spreading  angle. 

8.  One  can  model  the  governing  mechanism  in  these  flows  as  a  vortex-interaction 
process. 

9.  Acoustic  excitations  can  be  used  to  force  bifurcation  or  blooming  in  air  jets  at 
Reynolds  numbers  up  to  at  least  100,000  and  Mach  numbers  up  to  0.22.  With 
sufficient  power  in  the  forcing,  the  phenomena  can  probably  be  effected  at  any 
Reynolds  number  and  subsonic  Mach  number. 

4.  Future  plans 

Additional  experiments  in  the  high-speed  air  facility  are  needed  to  establish  the 
excitation  amplitudes  required  for  bifurcation  and  blooming  flow  control  of  high  Reynolds 
number  jets  at  moderate  Mach  Numbers.  This  is  now  being  accomplished  under  other 
AFOSR  funding. 


Methods  for  use  of  the  flow  power  to  produce  the  required  oscillation  need  to  be 
developed  if  the  phenomena  are  to  be  produced  in  large  jets  at  high  speeds.  Methods  for 
accomplishing  this  are  now  being  explored  under  other  AFOSR  funding. 
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abstract 

Thin  paper  la  concarnad  with  tha  control  of  a 
round  Jae  discharging  Into  quiescent  fluid  •  Tha  fluid 
la  water,  and  tha  Jae  la  azeltad  by  a  combination  of 
axial  and  orbital  axeleaelona.  Tha  axial  axcieaelona 
ara  lncroducad  through  notion  of  a  diaphragm  on  tha 
platon  that  drive*  tha  naln  flow,  and  tha  orbital 
axcltaelona  ara  producad  by  orbital  notion  of  tha  tip 
of  eh*  oozzla.  Tha  Say no Ida  o unbar  baaad  on  dlaaaear 
and  axlt  velocity  la  4300.  Sona  aurprlalng  rasults 
ara  obtalnad  whan  tha  ratio  (X)  of  tha  axial  frequency 
to  tha  orbl tal  frequency  la  la  tha  range  of  1-4 .  Whan 
R  "  Z,  tha  Jet  blfurcataa  Into  two  aeperate  ring- 
vortex  train*,  and  whan  R  -  3,  the  Jae  divides  into 
three  Jet*  In  an  equilateral  array  •  When  R  la  non- 
integer  In  tha  range  of  1.6-3 .2,  the  Jae  "bloom* '  and 
fill*  a  con*  of  80*  included  angle  with  vertex  ring* . 
Thla  produce*  a  vary  rapid  apraadlng  of  the  Jae  .  Tel¬ 
ocity  aaaaurananca  mad*  at  6 .3  and  8  dlanatara  away 
fron  the  noaal*  show  a  double-peak  profile  for  the 
bifurcating  and  blooning  Jata.  Thane  raaulta  ahow 
that  large  change*  can  ba  producad  In  round  Jet*  under 
controlled  excitation. 

SOMEWCLATTOE 

Orbital  axeltaelon  amplitude .  Peak  to  peak 
displacement  of  tha  noaal*  lip. 

OlaaMtar  of  eh*  noaal* . 

Axial  excitation  frequency,  Ha. 

Orbital  excltaelon  frequency,  a*. 

Ratio  of  axlal-eo-orbltal  excitation  frequency, 

W 

Reynold*  number  baaed  on  the  average  Jet  veloc- 
1*7  and  Jot  dlaneear,  °aT,t>/v. 


r  Radial  axla  away  from  tha  Jet  centerline  In  a 

cylindrical  coordinate  system. 

St(  Axial  Serouhal  number  baaad  on  diameter  and 
average  Jet  velocity,  faD/TJ4T- . 

St0  Orbital  Strouhal  number  baaad  on  dl am* tar  and 

average  Jat  velocity,  taDfxtt9m’ 

D  Mean  axial  velocity. 

°ava  Average  Jet  velocity  baaed  on  Jat  diameter  and 
the  volume  discharge  rata  from  tha  noaala . 

0Q  Mean  centerline  velocity  of  cha  natural  Jet  mea¬ 
sured  at  0.13  diameters  away  from  tha  noaala. 

u'  Soot -mean- square  value  of  tha  axial  velocity. 

x  Straamwlsa  direction  along  tha  cantarllo*  of  the 

Jat  la  a  cylindrical  coordinate  system. 

9  Aalnuthal  angle . 

v  Kinematic  vlacoalty. 

UTTRODCTCTIOH 

Numerous  studla*  In  tha  peat  have  triad  to  con¬ 
trol  tha  evolution  of  round,  turbulent  Jata  ualng 
controlled  excitation.  Sarohia  and  Baraal  [11  attend 
tad  to  Improve  mixing  by  pulsing  the  naan  flow  of  a 
round  turbulent  Jae.  Tlaea  [2]  advocated  the  us*  of 
oscillating  Jaes  end  developed  a  fluidic  aceuacor . 
Parry  and  Urn  [3|  obtained  long,  coherent,  vortical 
structures  while  exciting  a  co flowing  Jat  with  lateral 
oaclllatlona .  The  present  study  waa  Inspired  by 
Binder  and  Parra-Marinat  [6] ,  who  showed  that  tha 
effects  of  flapping  a  two- dimensional  Jae  are  felt  far  . 
away  from  tha  noaala .  Thla  lad  to  the  Idea  of  con-  j 
.  trolling  eh*  avolueloa  of  Jat*  through  combined  axial 
and  orbital  excltaelon. 

Thar *f ora,  tha  objective  of  thla  Investigation 
waa  to  study  tha  near- field  response  of  a  round  Jat 
under  varying  combinations  of  axial  and  orbital  exci¬ 
tation  frequencies.  Tha  raaulta  war*  documented  with 
eh*  aid  of  flew  visualisation  using  dye  Injection. 


1.7 


0.2 


Velocity  profile*  at  various  dowuacraaa  location*  wart 
aeasurad  with  a  Lasar-Dopplar  aneaonetry  syetan  haring 
aa  optieal-flbar  link  baevaan  the  laaar  and  the  crana- 
nlttlng  optica . 

tmtmwcAL  f  acuity 

Tha  facility  eoaalaead  of  a  large  plaeon-and- 
cy  Under  arrange*  sot  aa  ehowu  In  Fig.  1.  Tha  Jac  waa 
foraa J  by  pushing  wntar  oat  of  a  1/2-lnch-dlanacar  CD) 
nozzle  which  a Kl tad  Into  cha  hottan  of  a  t quart  plaxl- 
glas  cask.  Tha  placon  was  activated  aftar  cha  water 
waa  at  raat  and  fraa  of  any  notion .  Sinusoidal  par- 
curbaclona  vara  lapoaad  an  cha  aaan  flow  chrouth 
notion  of  a  flaxlbla  diaphragm  alcuacad  on  cop  of  tha 
platoa  •  Tha  dlaphraga  was  aac  Into  notion  by  a  hyd¬ 
raulic  ayaeaa  aad  driven  by  a  aocor  through  a  aat  of 
can pound  levers  aad  eaaa  •  Tha  dacalla  of  chla  ayacaa 
haw*  boon  deacrlbad  by  Bouchard  [3]. 
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■aaaura  the  atraaawlaa  eoapooanc  of  velocity.  Craan 
laaar  light  waa  guldad  to  tha  transmitting  opclcs 
(TSI,  Inc.)  through  a  polarlzatlon-praaarrlng,  slnglc- 
aoda  optical  fiber  (NIC  nodal  F-SPV-10)  a*  ahown  in 
Fig.  3.  Tha  optical  flbar  allowed  flexibility  In 
allganaat  aad  traweraa  of  cha  naaaurlng  voluaa .  The 
«f f active  naaaurlng  voluaa  waa  0J  n  wld*  and  2  no 
long.  Although  aany  ayateaa  have  bean  dealgned  and 
propoaad  baaed  on  optical  flbara,  cha  ayacaa  config¬ 
uration  follow*  that  of  Meakln  at  al .  [6]  for  alv- 
pllcley.  The  dacalla  of  chla  ayacaa  and  cha  orbital 
excitation  aachanlan  have  bean  deacrlbad  pravloualy 
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EXPERIfCNTAL  APPARATUS 


Fig.  1. 


tic  of  axperlanncal  facility  • 


Tha  orbital  axcltaelon  waa  achieved  by  novlng  cha 
tip  of  cb*  not* la  In  an  orbital  faahlon  through  a 
drlvn  plata  aad  four  equally  spaced  push-rods.  Aa 
shown  In  Fig.  2,  cha  peak— to— peak  notion  (Ag)  of 
cha  dp  of  cha  oossla  waa  approxiaacaly  0.918  Inches. 
The  spring-loaded  nosxle  chua  pracasaad  about  eha 
ayacaa  axis  at  an  angla  of  0.2*. 
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ESFElDOnrrAL  METHOD 


To  Tlaoallx*  tha  structure  la  tha  flow,  diluted 
food  coloring  waa  Injected  elrcuaferaatlally  through 
cha  dye  alota  shown  in  Fig.  2.  The  dye  Injected  Into 
tha  boundary  layer  of  the  nozzle  aarka  cha  shear  layar 
as  It  developed .  Both  33  m  anapahoca  and  18  an  clna 
flla  war*  aad*  of  tha  flow*. 

A  alngla-eonponant,  lasar-Dopplar  aoanonacry 
ayatan  In  tha  forvard-acatcar  aoda  waa  an  ploy ad  to 


Fig.  3.  Schanatle  of  lasar-Dopplar  ananoaacry  ayacaa. 
NATURAL  JET 

Ac  a  Reynolds  nuabar  of  *300,  cha  obaarvad  near 
flald  aeruceura  of  eha  natural  Jac  without  artificial 
axcltaelon  Is  eonslacanc  with  tha  obssrvatlona  of 
Browand  and  laufar  (8)  and  Bouchard  (3).  A  snapshot 
of  eha  natural  jat  using  dy*  injaetlon  can  be  found  In 


EFFECTS  OF  COMBINED  AXIAL  AND  ORBITAL  EXCITATION 


7i,  i.  Close  to  the  nozz  It ,  che  shear  layer  sepa¬ 
ratee  from  the  lip  and  rails  up  Into  vorcex  rings. 
The  formation  of  the  rings  Is  Irregular  In  else.  In 

present  flow,  the  average  shedding  race  was  26  Hz , 
haaad  oa  16  as  movies  .  The  ring  vortices  underwent 
•nceractlon  In  which  the  downstream  ring  would  grow 
lad  slow  down  while  the  following  ring  Increased  In 
speed  as  It  shrank  In  size.  The  upstream  ring  then 
passed  through  the  cencar  of  che  ring  In  front  and 
finally  coalesced  together  to  form  a  single  ring.  At 
ibouc  five  diameters,  the  rings  broke  down  into  chree- 


7ig  .  ■>.  Natural  jet  at  Re  •  4300. 


Telocity  measurements  of  the  nacural  jet  taken  at 
3.15  diameters  downstream  fx/D  •  0.15)  showed  chat 
the  aean  velocity  was  flat  within  1.51  (Fig.  5).  The 
average  jet  exit  velocity  H,,,)  «*•  0-207  m/s,  and 
the  centerline  velocity  !"  )  was  0.327  m/s.  The 
shear  layer  had  a  half-slope  thickness  of  0.07  diam¬ 
eters  with  a  centerline  fluctuation  Cu')  of  less  chan 
'.I  at  X/D  *  3.15.  Telocity  profiles  of  the  nacural 
jet  jo  to  eight  diamecers  downstream  appear  symmetric, 
with  s  smooth ,  bell-shaped  curve . 

1  .'c" 
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rig.  3.  >Loc;ty  profllas  cf  natural  Jet. 


In  the  present  study,  a  sinusoidal  pertubacioo  of 
171  was  Imposed  on  che  mean  flow.  The  axial  Jac 
responded  to  che  Imposed  excitation  by  shedding  vortex 
rings  at  che  excitation  frequency.  Visualization 
showed  chat  all  vortex-ring  coalascanca  was  suppressed 
by  this  excitation. 

Uhen  the  axial  and  che  orblcal  excitations  were 
borh  Imposed  on  che  jet,  dramatic  changes  were  ob¬ 
served.  The  dimensionless  number  R  defined  as  the 
ratio  of  axlal-co-orbltal  exclcacloo  fraouancy  la  an 
Important  paraaacer  characterizing  che  flows.  An 
interesting  overview  of  the  various  cypes  of  Jscs  can 
be  obealned  when  the  axial  excltadon  freauency  (f  ) 
la  fixed  at  12  Hz  while  increasing  the  orbital  excita¬ 
tion  freauency  (f0)  as  R  decreases,  where  Re  ■ 
4300.  This  will  now  be  reviewed. 

At  low  orblcal  frsausncles,  the  global  features 
of  che  flow  are  similar  to  chose  of  an  axially  excited 
jet.  Vortex  rings  are  observed  co  shed  regularly  at 
the  axial  excitation  frequency  of  12  Hz.  In  addition, 
the  vortex-ring  train  exhibits  a  slightly  sinuous 
behavior. 

When  che  orbital  frequency  is  gradually  in- 
creaaad,  the  vortax  rlnga  begin  co  tilt  and  collide 
with  each  other,  breaking  up  at  about  five  diameters 
from  the  nozzle.  During  this  process,  tha  dys  which 
marks  tha  shsar  Isvar  is  obsarvad  to  disperse  ae  a 
larger  angle  when  enmnared  with  a  natural  Jat.  The 
obsarvad  dispersion  angle  la  about  40*  whan  R  •  4. 
When  R  is  further  decreased  co  a  value  between  3  and 
4,  vortex  rings  lntermiccendy  move  away  from  che  jet 
axis  at  an  angle.  The  vortex  rings  that  move  away 
from  cha^  jat  axis  remain  coherent  for  five  co  ten 
diamecers. 

When  che  value  of  R  is  aooroximacsly  3,  the  jac 
splits  up  into  three  branches  with  a  slight  preces¬ 
sion.  The  orecession  disappears  when  R  -  3,  with 
the  jac  forming  an  aqullataral  array.  Because  of  this 
flow  pattern,  che  flow  is  called  a  ’trifurcating'  jet. 
A  similar  phenomenon  can  be  observed  when  R  i  2 
(Fig.  6).  In  this  case,  Che  Jet  hlfurcates  tnco  s  T- 
shaped  pattern  in  opposite  directions.  Flgura  '  snows 
an  on-axis  view  of  tha  ’bifurcating’  Jat  issuing  uuc 
of  tha  plana  where  R  *  1.9.  The  precession  of  the 
bifurcating  Jac  la  tndlcaced  by  evo  spiraling  trails 
of  dys  deposited  by  che  vortex  rings  as  chey  breax 
up.  Occasionally,  a  secondary  forking  of  one  of  che 
branches  can  be  seen.  This  type  of  Jec  is  called  a 
’bifurcadng’  Jet  because  of  che  unique  forking  beha¬ 
vior  of  tha  vorcsx  rings.  The  angle  between  che  axes 
of  che  two  branch  of  vortex  rings  nas  oeen  observed  co 
be  a a  large  as  60*. 

Ac  R  “  1.7,  vortax  rings  can  be  seen  living 
away  from  che  centerline  along  a  side  of  a  cone  In  all 
directions.  A  plccura  of  this  phenomenon  is  shown  in 
Fig.  3.  The  flow  field  resembles  chac  of  a  flower 
blooming  when  viewed  on  axis  (Fig.  9),  and  hence  che 
flow  Is  called  a  ’blooming’  jec.  This  phenomenon 
persists  for  a  range  of  R  *  3. 2-1. 6. 

When  the  orblcal  exclcadon  freauency  Is  further 
Increased  such  chac  R  falls  below  1.6,  che  disper¬ 
sion  angle  of  che  vortex  rings  drastically  decreases. 
At  R  <  1.3,  che  vortex  rlnga  break  up  vary  earlv 
during  the  roll-up. 
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CLASSIFICATION  of  flow  field 
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Fig.  *.  End  view  of  bifurcating  Jet. 


Fig.  A.  Blooming  <ec  at  R  *  1 .7  . 


Fig.  a.  Zri  view  o.‘  hlooalae  jet  • 


In  an  attcape  to  map  out  tha  diffarant  types  of 
Jeca,  dye  lnjacclon  uaa  uaad  co  visualize  various 
flows  at  Reynold!  numbers  of  2800,  4000,  3500,  8300 
and  10,000.  Soma  of  tha  Jett  are  marked  on  tha  paraa- 
ecer  plot  shown  in  Fig.  LO.  Tha  frequencies  are  nor- 
aalizad  by  D  and  .  0  to  fora  axial  (Sta)  and 
orbital  (St0)  Strouhal  nuaban .  Tha  open  jynbols 
repreaent  excited  jet*  whose  vortex  rings  have  been 
observed  to  break  up  before  being  able  co  propagate 
away  froa  the  centerline .  The  closed  symbols  repre¬ 
sent  exclced  jets  whose  vortex  rings  remain  coherenc 
acre  chan  two  diameters  away  from  the  centerline. 
Within  the  region  roughly  marked  by  the  closed  sym¬ 
bols,  bifurcating,  trifurcating,  and  blooming  jets  can 
be  found.  For  simplicity,  the  area  is  called  the 
‘blooming*  region.  The  region  lies  roughly  between  R 
•  1.6-3 .2.  Furthermore,  the  region  is  confined  be¬ 
tween  St,.  “  0.33-0.75  and  Sc  -  0.15-0.35.  When  S 
-  2,  the  jet  bifurcated  Into  a  T-shaped  pattern,  and 
when  R  •  3,  the  jet  aplita  up  Into  three  branches. 
If  the  value  of  I  !  1  or  3,  the  jet  la  observed  to 
preceaa  with  the  corresponding  flow  patterns.  When 
R  takes  on  non-lntagar  values  within  the  region,  the 
Jets  “bloom*  to  fill  a  cone  of  30*  Included  angle. 


Fig.  10.  Parameter  plot  of  bifurcating  and  blooming 
Jets . 

FLOW  STSOCTTOES 

Baaed  on  the  dye-lnjectlon  photographs,  the  ap¬ 
pearance  of  che  excited  jets  la  very  different  from 
that  of  the  natural  jet .  Besides  che  tremendous 
lncreaea  in  the  spreading  angle,  a  bifurcating  Jec 
loses  its  axisymmecry  and  exhibits  two  planes  of 
symmetry.  Tha  first  plane  of  symmetry  lias  along  a 
diametric  plane  containing  the  two  branches  of  the  jec 
called  the  plane  of  bifurcation.  The  second  plane  of 
symmetry  bisects  Che  two  branches  and  is  called  the 
bisecting  plane.  La  the  case  of  the  blooming  jet,  the 
axisymmecry  appears  co  be  racalned. 

Within  these  complicated  flow  fields,  the  follow¬ 
ing  fluid  motion  can  be  deduced  .  Figure  11  shows  a 
schematic  of  a  bifurcating  Jtt  in  which  the  solid  Xlna 
marks  tha  boundary  between  the  jec  and  ambient  fluid  . 
In  che  first  few  diameters,  layers  of  jec  and  aablenc 
fluid  roll  up  to  form  vortex  rings.  Dosmstrsaa  of  the 
fork  In  the  flow,  the  jet  fluid  must  be  stretched  back 
and  forth  between  two  branches  because  of  che  alter¬ 
nating  motion  of  suecaaaive  rings.  In  the  same  fash¬ 
ion,  blooming  jets  would  exhibit  s  similar  moeloo  in 
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*hie>i  eh*  Jet  fluid  1*  stretched  back  and  forth  in  *11 
ssiauthel  direction*. 


ns-  11.  Schematic  of  *  bifurcating  Jet  with  trajec¬ 
tories  of  Indue ad  downwaah. 

tklocitt  profits 

Mean-velocity  prof  Hat  measured  at  various  down- 
•eream  location*  of  cha  blooalng  jac  (1  •  2.4)  and 
la  the  plana  of  bifurcation  of  tha  bifurcating  Jet  (R 
•2)  rataa I  some  lncaraaclng  developments .  Tha  Ray- 
oolda  auabar  t*  4300  and  f4  ■  12  Sz .  In  tha  first 
chraa  dlaaacsrs,  both  bifurcating  and  blooming  Jaca 
exhibit  a  flat  coyhae-ahapad  proflla  which  gradually 
•mooches  out  Into  a  ball-ahapad  curve  Ilka  a  natural 
Jec .  from  four  to  eight  diameters  downstream,  a 
double-peek  proflla  data  lope  In  both  bifurcating  (Fig. 
12)  and  blooming  (Fig.  13)  Jets.  Thai  a  peaks  appear 
to  coincide  with  tha  location  of  tha  coherent  tor tea 
ring*  which  split  up  and  more  away  from  the  Jet 
centerline.  Since  tha  vortex  rings  ara  probably  tha 
aoet  energetic  strueturea  la  these  flows,  this  obser- 
tation  ts  ex  pact  ad .  what  la  surprising  are  tha  addi¬ 
tional  dips  within  tha  double-peak  profiles  of  tha 
bifurcating  jet  et  6.3  and  3  disasters.  Tha  double 
peeks  of  rha  blooming  Jet  develop  at  about  tha  same 
location  but  without  tha  additional  dtps. 

In  an  attempt  to  understand  cha  origin  of  cha 
dlaa,  valoelty-tism  trscaa  for  one  second  at  x/D  ■ 
6  o  ware  recorded  on  in  oecilloscope  at  locationa 
Inside  tha  dip  and  *e  tha  peaka  of  tha  bifurcating 
Jet .  Figure  14  show*  a  comparison  between  the  two 
cr*caa>  Although  both  traces  reveal  a  regular  fluc¬ 
tuation  of  6  %  (vortax-paaaing  frequency)  and  about 
“a*  peak  velocity,  tha  boetom  trace  recorded  In 
the  dip  shows  such  lower  valleys  .  If  tha  schematic  of 
a  bifurcating  Jet  is  analyzed,  as  shown  in  Figure  ll, 
four  trajectories  can  be  seen  along  which  a  large 
downwaah  component  of  velocity  la  induead  by  cbe  rr* 
passing  of  vortex  rings .  Tho  repeated  psaaago 
bf  the  vortex  rings  in  this  location  could  load  to  tha 
formation  of  the  dip#  in  the  velocity  profiles  as 
observed.  In  the  blooming  Jet,  dlpe  do  not  appear 
because  the  ringe  do  not  follow  each  other  along  flxad 


crajeecorlea  as  in  a  bifurcating  Jet.  Therefore,  the 
downwaah  component  of  velocity  induced  by  aach  vortex 
ring  la  not  reflected  in  the  time-averaged  velocity 
profiles . 
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-Fig.  12.  Mean-velocity  profiles  of  bifurcating  Jet. 
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Fig.  13.  Mean-velocicy  profile*  of  blooming  Jet. 


Velocity  Trace  of  Bifurcating  Jet(Qxtf) 
at  X/D  •  6.5  for  1  sec. 
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Fig  .  14  .  Velocity  tracel  of  a  bifurcating  jet . 
CONCLUSIONS 

In  ausTsary,  the  spreading  angle  of  a  round  jet 
haa  been  tremendously  Increased  by  imposing  axial  and 
orbital  excitation  on  the  flow.  The  axial  excitations 
are  small  pulsations  of  the  seen  flow,  and  the  orbital 
excitations  are  small  orbital  motions  of  the  Up  of 
the  oorile.  Different  flow  flelde  called  bifurcating 

and  blooming  Jecs  have  been  observed .  A  bifurcating 

Jet  conelsca  of  vortex  rings  split  up  into  two  trains 

In  a  Y-shaped  natters  and  la  produced  lAien  Che  axial 
excitation  frequency  la  two  times  that  of  Che  orbi¬ 
tal.  When  the  ratio  of  excitation  frequencies  takas 
an  aoo-taceger  values,  the  Jet  'blooms*  to  fill  a  cone 
of  80*  Included  angle  with  vortax  ring*  moving  along 
che  aide.  Obeervaclon*  based  on  dye  Injection  and 
velocity  measurements  show  thee  controlled  axial  and 
orbital  excitation  can  produce  a  large  change  in  che 
development  of  round,  turbulent  Jecs. 
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— Abstract - 

Three-dimensional  vortex  simulations  are  used  to  study  the  remarkable 
structure  of  the  round  jet  with  imposed  axial  and  orbital  excitations  and 
to  investigate  the  dynamics  of  an  inhomogeneous  mixing  layer  in  which 
the  unperturbed  state  is  an  infinite  array  of  Stuart  vortices.  The  simula¬ 
tions  reveal  the  nonlinear  vortex  interactions  of  importance  that  control 
the  structure  and  development  of  these  flows. 
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Abstract 

Three-dimensional  vortex  simulations  are  used  to  study  the  remarkable  structure  of  the 
round  jet  with  imposed  axial  and  orbital  excitations  and  to  investigate  the  dynamics  of  an 
inhomogeneous  mixing  layer  in  which  the  unperturbed  state  is  an  infinite  array  of  Stuart 
vortices.  The  simulations  reveal  the  importance  of  the  nonlinear  vortex  interactions  that 
control  the  structure  and  development  of  these  flows. 

1.  INTRODUCTION 

In  this  paper  we  discuss  the  computer  simulation  of  complex  three-dimensional  fluid  flows 
by  vortex  methods  and  give  two  applications.  By  simulation  we  mean  that  a  large  range  of 
the  important  scales  of  motion  in  space  and  time  are  computed  directly,  as  in  large-eddy 
simulations.  The  flows  we  consider  are  at  constant  density  and  have  high  Reynolds 
numbers. 

In  a  three-dimensional  vortex  method,  vector  elements  of  vorticity  move  at  or  near  the  local 
fluid  velocity  with  the  vectors  strained  by  the  local  velocity  gradient.  These  elements  are 
required  only  where  the  vorticity  is  nonzero  and  are  created  only  (1)  to  represent  initial 
fields  of  vorticity  or  vorticity  injected  through  computational  boundaries.  (2)  to  satisfy  no¬ 
slip  boundary  conditions  or  (3)  during  "grid"  refinement  when  one  element  may  be  split  into 
two  or  more. 

In  many  applications,  the  vorticity  field  is  well  represented  by  isolated,  thin  tubes.  If  the 
tubes  are  assumed  to  have  zero  cross  section,  a  logarithmic  divergence  in  the  induced 
velocity  is  encountered.  Thus  the  finite  size  of  the  vortex  core  must  be  taken  into  account 
in  developing  equations  for  vortex  tube  dynamics.  As  a  result  of  such  analyses,  accurate 
equations  for  the  long-wavelength  dynamics  of  thin  vortex  tubes  are  available  and  form  the 
basis  for  a  three-dimensional  vortex  method,  as  discussed  in  Section  2.  To  discretize  the 
vortex  tubes,  the  space  curve  of  each  tube  or  filament  is  marked  with  a  sequence  of  node 
points.  Velocities  of  all  the  N  nodes  on  the  vortex  fiiameiits  are  computed  from  the  infinite 


medium  Green's  rV~tion  (Biot-Savart  law),  requiring  0<~XZ)  '■'erauons  per  time  step.  The 
strain  on  the  vortex  elements,  which,  roughly  speaking,  are  thw  vectors  connecting  adjacent 
nodes  points  on  each  filament,  is  accounted  for  automatically  by  the  method. 

In  Section  3  we  discuss  the  simulation  of  the  excited  round  jet.  using  the  thin-filament 
approximation.  In  this  application,  it  had  been  observed  experimentally  that  dramatic 
changes  in  the  structure  of  the  jet  fluid  are  achieved  by  introducing  axial  and  orbital 
excitations  [61.  The  interplay  between  computation  and  ongoing  laboratory  experiments  has 
led  to  a  basic  understanding  of  the  mechanisms  that  govern  the  structure  of  the  jet  and  how 
they  might  be  controlled. 


In  other  applications,  the  vorticity  field  is  a  relatively  smooth  continuum.  In  this  situation, 
vortex  elements  must  be  densely  packed  to  represent  properly  the  dynamics  of  the 
continuum.  Recent  mathematical  proofs  have  shown  that,  when  properly  constructed, 
vortex  methods  will  produce  solutions  that  converge  to  the  solution  of  the  Euler  equations 
as  the  number  of  computational  elements  increases  [1,2,5].  However,  because  a  relatively 
large  number  of  elements  is  required  for  these  applications,  the  cost  of  computation  can  be 
very  high  if  the  dynamics  are  computed  via  the  Biot-Savart  law.  Alternatively,  the  vortex- 
in-cell  technique,  described  in  Section  4,  may  be  used  to  compute  the  vortex  interactions 
thereby  significantly  reducing  the  price  of  the  simulation.  In  Section  5  we  use  the  vortex- 
m-ceil  method  to  study  the  dynamics  of  perturbed  Stuart  vortices  -  a  model  of  the  mixing 
layer  witn  two-dimensional  structure. 

2.  DYNAMICAL  EQUATIONS  FOR  VORTEX  FILAMENTS 

For  a  constant  density  flow,  the  dynamics  of  the  vorticity  field  at  is  given  by 

-§2L  3  +  *  *-Vu  +  vVzat  ,  (1) 

Dt  St 

where  u  is  the  velocity  and  v  is  the  kinematic  viscosity.  From  the  relations  V  u*0  and 
V  x  u  —  at  we  find  that 

V2u  -  -7x«  ,  (2) 


so  that  u  may  be  expressed  as 


u(r,t) 


1  f  (r  -  rO  x  M(r',t) 

4fl-J  ir-r'!3 


dr*  ■)•  V  , 


(3) 


where  r  is  the  space  vector  (x,y,z)  and  d>  is  the  potential  associated  with  the  homogeneous 
solution  of  12)  required  to  satisfy  boundary  conditions. 


For  numerical  purposes  we  assume  that  the  vorticity  may  be  represented  as  a  system  of 
vortex  filaments  following  the  three-dimensional  space  curves  rj(f,t),  (i  *  1,  2,  3,  ...) 
where  f  is  a  parameter  along  a  given  curve.  The  velocity  field  is  then  given  by 


(4) 
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r,  C(r-fi(f',t))  x  (8r,/8f)g(|r-r,|/(7,) 
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where  g  is  a  smoothing  function  that  depends  on  the  assumed  structure  of  the  vortex  core 
as  discussed  in  (8],  For  the  study  described  in  the  next  section  we  used  the  function 


g(c) 


c3 

(c ‘  4-  a)3/2 


(5) 


in  (4)  with  a  — 0.413,  corresponding  to  Gaussian  cores  [7],  and  computed  the  velocities  of 
the  node  points  at  fk,  (k  —  1,  2,  3,  ...)  using 


3r(fk,t) 

3T“ 


u(r(fk,t))  . 


(.6) 


Viscous  effects  are  neglected  except  for  the  creation  of  vorticity  at  a  no-slip  boundary. 


A  simple  way  to  approximate  the  geometry  of  each  space  curve  is  to  use  a  sequence  of 
linear  segments  between  successive  nodes.  With  this  approximation  we  can  obtain  an  exact 
analytical  expression  for  the  integral  in  (4)  with  g  given  by  (5).  In  an  improved  method 
that  we  recently  developed,  the  trapezoid  rule  is  used  to  integrate  along  each  space  curve 
with  parametric  cubic  splines  to  estimate  the  required  derivatives  3r/3f  at  each  node.  The 
new  method  requires  only  half  the  number  of  arithmetic  operations  per  node  point  and, 
generally,  only  half  the  number  of  node  points  for  equivalent  accuracy.  As  an  illustration, 
the  convergence  of  the  self-induced  velocity  of  a  vortex  ring  are  shown  in  Figure  1  for  the 
linear-segment  scheme  and  the  spline-based  method  as  a  function  of  the  number  of  points. 
As  a  consequence  of  the  savings  in  operation  count  and  number  of  node  points  (for  the 
same  accuracy),  the  spline-based  scheme  is  approximately  a  factor  of  eight  faster. 

3.  SIMULATION  OF  THE  PERTURBED  ROUND  JET 

In  a  typical  unperturbed,  turbulent  round  jet,  the  shear  layer  exiting  the  jet  nozzle  rolls  up 
into  distinct  vortex  ring  structures  that  eventually  merge  and  appear  to  loose  their 
coherence.  In  the  experimental  work  of  Lee  &  Reynolds  [6]  it  was  observed  that  small 
perturbations  can  significantly  alter  the  development  of  round  turbulent  jets.  They 
introduced  both  axial  and  orbital  excitations.  By  pulsing  the  jet  axially,  they  control  the 
frequency  at  which  the  vortex  ring  structures  are  formed.  The  orbital  excitation  is  produced 
by  moving  the  nozzle  tip  in  a  circular  orbit  about  the  nominal  centerline  of  the  jet.  If  the 
orbital  period  rfa  small  integer  (*  multiple  of  the  period  of  the  vortex  ring  formation, 
the  jet  splits  into  nr  distinct  trains  of  vortex  rings  equally  spaced  around  the  azimuth.  For 
Or  — 2.  the  angle  between  the  branches,  referred  to  as  the  bifurcation  angle,  was  found  to 
vary  from  40°  to  80®,  depending  on  the  Sirouhal  number  (St  »  fD/<U>). 


In  this  section  we  investigate  the  dynamics  of  this  process  for  nr“2,  the  bifurcating  jet, 
using  a  three-dimensionai  vortex  simulation  with  the  thin-filament  approximation.  The 
important  vorticity  in  the  flow  is  assumed  to  be  in  the,  form  of  a  semi-infinite  cylindrical 


sheet  of  azimuthal  vorucity,  representing  the  source  flow,  and  a  collection  of  vortex 
filaments  as  illusti  d  in  Figure  2.  The  velocity  field  mdu.  j  by  the  vortex  sheet  was 
defined  analytically  and  an  efficient  scheme  was  developed  for  its  numerical  evaluation. 
New  vortex  rings  are  created  periodically  at  the  jet  exit  with  a  circulation  consistent  with  the 
flux  of  circulation  given  by  dT/dt»  0.5U-*  where  L’c  :s  the  jet  centerline  velocity.  The 
effect  of  orbital  excitation  is  approximated  in  the  simulation  by  rotating  the  cylindrical 
vortex  sheet  in  a  circular  orbit  about  the  nominal  centerline  of  the  flow. 

Typical  experimental  and  computational  results  are  shown  in  Figures  3  and  4  for  a  case  of 
large  bifurcation  angle.  In  the  experiment,  dye  was  injected  circumferentially  at  the  jet  exit, 
to  marie  the  vortex  structures  in  the  jet  fluid.  The  computational  results  have  shown  that 
the  axial  spacing  between  rings,  which  depends  inversely  on  St,  is  a  crucial  parameter.  For 
the  spacing  that  yields  maximum  spreading,  the  successive  rings  strongly  interact  at  a  certain 
downstream  location,  launching  one  another  along  widely  separated  paths  (Figure  4b).  If 
the  spacing  is  too  large,  the  interaction  is  too  weak  to  produce  wide  angles  (Figure  4a).  If 
the  spacing  is  too  small,  successive  rings  collide  and  merge  (Figure  4c).  As  a  consequence, 
one  obtains  a  bifurcation  angle  that  increases  with  St,  but  with  low  and  high  frequency 
cutoffs  as  shown  in  Figure  5a.  The  experimental  results  [6]  presented  in  Figure  5b  show  a 
similar  relationship  between  bifurcation  angle  and  Strouhai  number.  However  the 
maximum  angles  achieved  in  the  experiment  are  twice  those  of  the  simulation,  and  the 
excitation  amplitude  is  an  order  of  magnitude  less.  We  are  currently  exploring  a  number  of 
possible  reasons  for  these  discrepancies.  Errors  due  to  the  thin-filament  approximation  may 
be  significant  and  the  details  of  the  flow  near  the  nozzle  may  need  more  accurate 
representation.  Nevertheless  the  vortex  simulations  predict  the  correct  trends  and  we 
beiieve  that  they  demonstrate  that  vortex  interaction  is  the  key  mechanism  governing  the 
structure  of  the  flow. 

4.  VORTEX-IN-CELL  METHOD 

As  the  number  of  vortex  node  points  increases  beyond  a  few  thousand,  the  vortex  method 

described  in  Section  2,  a  direct-interaction  scheme  using  the  Biot-Savart  law.  becomes 

prohibitively  expensive,  even  on  presently  available  supercomputers.  In  this  section  we 

discuss  an  alternative  scheme  to  compute  node  velocities  at  a  much  lower  cost.  The  idea  is 

to  retain  the  Lagrangian  treatment  of  the  vorticity  field  but  to  solve  the  Poisson  equation  for 

the  velocity  field  on  a  fixed  Eulerian  mesh.  By  use  of  fast  Poisson  solvers  on  a  mesh  of  M 

grid  points,  the  operation  count  for  this  step  can  approach  O(MiogM).  Additional  required 

steps  are  (l)  the  generation  of  mesh  values  of  voracity  from  the  Lagrangian  represenution 

and  (2)  interpolating  velocities  from  the  mesh  back  onto  the  Lagrangian  points.  Therefore 

the  computing  time  per  time  step  for  the  vortex-m-cell  method  is  O(N)  i-  O(MlogM).  a 

substantial  reduction  over  the  0<N:)  time  required  for  the  direct-interaction  scheme. 

However  the  vortex-in-cell  scheme  requires  a  more  complex  algorithm  and  the  introduction 

of  another  layer  of  numerical  approximation.  In  addition,  some  boundary  conditions,  such 

as  the  infinite  domain,  are  more  difficult  to  implement,  while  the  periodic  boundary 

« 

condition,  for  example,  is  simpler. 


In  the  three-dimen«tonai  vortex-m-cell  method  of  Couet  et  f3] ,  developed  for  flow  fields 
that  are  periodic  ..  ail  three  directions,  the  vorticity  field  is  .^presented  as  a  collection  of 
filaments  with  Gaussian  cross  section. 


«(r.t) 
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A  specially  developed  quadratic-spline  interpolation  scheme  is  used  to  minimize  the  errors  in 
deposition  of  vorticity  onto  the  mesh  and  in  interpolation  of  velocities  onto  the  Lagrangian 
points.  The  method  was  used  in  the  computation  of  model  problems  and  the  results  were 
compared  with  solutions  icnown  to  be  accurate  [3].  Excellent  results  were  achieved.  In  the 
application  below  we  use  the  extension  of  the  method  to  the  infinite  domain  [4],  suitable  for 
slabs  of  vorticity,  having  finite  thickness  in  one  direction. 


5.  DYNAMICS  OF  AN  INHOMOGENEOUS  MIXING  LAYER 


We  apply  the  vortex-in-cell  method,  described  in  the  previous  section,  to  the  study  of  the 
dynamics  of  perturbed  Stuart  vortices  [12].  The  vorticity  field  for  the  unperturbed,  two- 
dimensional  flow  is  periodic  in  the  stream  (x)  direction  with  period  L  and  is  given  by 


a»j(x.y) 


—  2ir  (1  —  o2)  /L _ 

[cosh(2iry/L)  -p  cos(2irx/L)]2 


(8) 


where  p€  [0,1]  is  a  parameter.  If  p  —  0,  a  parallel  tanh  shear  profile  is  obtained  and  if  p—  1 
one  obtains  a  row  of  point  vortices.  At  y  —  m  oo  the  streawwise  velocity  is  2: 1 .  Thus  the 
circulation  of  each  vortex  is  2L.  In  the  following,  we  study  the  case  p  —  0.25  and  set  L  * 
16. 


To  simulate  the  dynamics  of  this  vorticity  continuum  we  use  an  array  of  16x25  vortex 
filaments  for  each  Stuart  vortex.  Thus  the  spacing  between  filaments  is  Ax  — Ay  —  1.0.  In 
the  span  (z)  direction.  32  node  points  are  used  to  define  the  space  curve  of  each  filament. 
To  allow  for  the  study  of  subharmonic  disturbances,  the  computational  domain  includes  two 
Stuart  vortices  or  300  filaments  and  a  total  of  25600  node  points  in  a  32x3  Eulerian  mesh. 

Each  vortex  filament  represents  a  vorticity  field  given  by  Eq.(7)  with  an  effective  core 
parameter  <r.  Therefore  the  circulation  of  each  filament  may  be  determined  to  minimize 
square  error  in  representing  the  exact  vorticity  field  given  by  (8).  Figure  6  shows  the 
quality  of  the  representation  that  is  obtained.  Contours  of  vorticity  field  given  by  (8)  are 
plotted  together  with  the  field  obtained  numerically;  they  are  virtually  indistinguishable. 
The  rms  error  normalized  to  <m{>>1  is  only  0.12%.  The  exact  unperturbed  state  is  of 
course  steady  but  the  approximate  vortex  filament  representation  produces  a  time- 
dependent  flow.  Figure  7  shows  the  evolution  of  the  unperturbed  state  after  380  time  steps. 
Nearly  neutral  oscillations  are  observed  with  normalized  rms  error  equal  to  29%.  However, 
the  time-averaged  flow  remains  close  to  the  exact  steady  state  <  Figure  8)  with  normalized 
rms  error  equal  to  5.6%. 


Next  we  consider  sr  three-dimensional  perturbations  to  thr  uart  vortices.  The  linear 
stability  calculations  ot  Pierrehumbert  &  Widnaii  [10|  ire  available  for  comparison.  In  our 
calculations  the  vortices  are  centered  at  x-0  and  x  —  16  <—  — 16).  Thus  by  displacing  the 
vortex  node  points  in  the  x  -  y  plane  at  t  »  0  by  the  amount 

Sx-Sy-.coJ^-  coJ^.  . 
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where  «  is  a  small  parameter,  we  obtain  the  pairing  instability  for  a  —  1  and  the  translative 
instability  fora^O,  as  studied  in  [10]  for  spanwise  wavenumbers  0  —  0.1.2,  ■  •  • 

We  expect  that  after  an  initial  period  of  adjustment  the  most  amplified  mode,  if  one  exists, 
will  dominate  the  perturbation  for  intermediate  times  and  yieid  an  exponential  growth  rate. 
At  later  times  nonlinear  effects  should  play  a  role  in  the  dynamics.  These  features  are 
displayed  in  Figure  9  where  we  show  the  energy  in  mode  /3  —  1  for  the  initial  disturbance 
«*■  1,  /3*«  1.  In  Figure  10  we  plot  the  exponential  growth  rates  for  this  case  and  several 
others  along  with  the  results  of  Pierrehumbert  &  Widnaii  [10].  The  agreement  for  small  0 
is  good  to  within  a  few  per  cent.  Naicamura  et  al  [9]  observed  errors  of  a  similar  magnitude 
in  their  study  of  the  growth  rates  of  0  ■*  0  perturbations  to  parallel  shear  flows  using  a  two- 
dimensional  vortex  method  with  similar  resolution.  In  the  translative  instability  (Figure 
10a  ),  the  agreement  for  larger  0  is  not  good  quantitatively  but  shows  the  same  leveling  off 
to  a  broad  plateau.  For  the  a  ™  1  mode  (Figure  10b)  Pierrehumbert  &  Widnaii  observed  a 
high  0  cutoff  to  the  pairing  instability,  whereas  the  present  results  suggest  a  broader  domain 
of  instability.  One  obvious  possibility  for  these  discrepancies  is  insufficient  resolution  in  the 
present  results  which  causes  a  twofold  problem.  First ,  the  unperturbed  state  does  not 
consist  of  the  exact,  steady  Stuart  vortices,  but  is  time-dependent  (recall  Figure  7)  and. 
second,  the  resolution  of  the  perturbation  state  must  rely  on  the  relatively  coarse  grid  that  is 
used.  Another  possibility  is  that,  in  our  simulation;  a  branch  to  another  family  of  unstable 
modes  has  taken  place  as  0  is  increased  -  a  family  that  has  more  complex  structure  in  the 
vortex  core.  Finite  area  vortices  (ellipses),  when  subjected  to  an  external  strain  field,  are 
subject  to  an  infinite  hierarchy  of  instabilities  as  /3—*oo  [11,131.  If  the  strain  is  low,  the 


unstable  domains  in  0  are  separated  by  neutrally  stable  domains.  For  higher  strain,  the 
unstable  domains  completely  overlap.  Further  work  is  needed  to  establish  whether  such  a 
hierarchy  exists  for  the  Stuart  vortices.  We  suspect  that  it  does  and  that  we  have 
encountered  such  a  higher  mode  instability,  but  at  the  possibly  incorrect  wavenumber.  In 
Figure  11  we  show  contours  of  a where  at ^  is  the  computed  vorticity  field  for  the 
unperturbed  initial  condition  (e-*0).  Figure  l la  illustrates  the  mode  shape  for  the  pairing 
instability  a**  1,  0  —  0  after  380  time  steps.  The  mode  profile  agrees  qualitatively  with  that 
given  in  [10]  for  the  same  spanwise  wavenumber.  Figure  lib  illustrates  the  mode  shape  for 
the  translative  instability  a  —  0,  1;  the  structure  here  appears  more  complex  but  the 

results  are  not  conclusive. 


Of  course  the  primary  advantage  of  the  present  method  ss  that  fully  nonlinear  studies  of 
three-dimensional  inviscid  flows  are  possiibte.  Figure  12  shows  the  energies  of  the  various 
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spanwise  modes  in"'  the  nonlinear  regime  for  the  cases  a  —  0  3  —  2  and  for  a  -  1,  0  —  1. 
As  seen  from  the  . iSuaiizauon  of  the  vortex  filaments  in  Figure  13,  the  vortex  cores  are 
starting  to  overlap  when  t— 300  and  the  vortex  shape  deviates  noticeably  from  a  sinusoid 
when  t  — 400.  Studies  of  more  complex  physical  problems  are  planned,  including  vortex 
ring  interactions  witn  Stuart  vortices. 

8.  CONCLUSIONS 

Three-dimensional  vortex  methods,  employing  Lagrangian  elements  of  vorticuy,  are  well 
suited  to  the  study  of  many  complex  flows  at  high  Reynolds  numbers.  When  the  flow 
consists  mainly  of  isolated  tubes  of  vorticity,  a  single  computational  filament  may  be 
assigned  to  each  physical  vortex  and  the  dynamics  computed  via  the  Biot-Savart  law  at  a 
reasonable  cost.  Simulations  of  the  perturbed  round  jet  were  performed  with  this  scheme 
and  have  shown  that  three-dimensional  vortex  interactions  between  successive  vortex  rings 
are  responsible  for  the  observed  bifurcation  of  the  jet  Quid. 

If  the  physical  vortex  tubes  undergo  collision,  merging,  or  any  other  process  in  which  the 
structure  of  the  vortex  core  is  important  to  the  dynamics,  a  large  number  of  computational 
elements  must  be  used  to  represent  each  vortex.  In  many  cases  the  vortex-in-cell  method 
may  be  used  to  reduce  substantially  the  cost  of  computation.  In  this  paper,  the  vortex-in- 
ceil  method  was  used  to  study  the  linear  and  nonlinear  three-dimensional  vortex  interactions 
of  perturbed  Stuart  vortices,  a  fiow  field  that  is  representative  of  an  inhomogeneous  mixing 
layer. 
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Figure  8  Contours  of  analytic  vorticity  and  their  2-D  projection  together  with  the  time- 
averaged  computed  vorticity  (dashed)  after  384  time  steps.  Increment  in  2-D 
contour  levels  *  0.05. 


n 

*7o 


Ed  2 


0.0  100.0  200.0  300.0  400.0  500.0  600.0 

TIME 


Figure  9  Time  evolution  of  total  energy  in  mode  & 1  lor  the  initial  disturbance 
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Figure  10  Growth  rate  versus  (3  tor  (a)  a  *  0  and  <b)  a**  1  Solid  line*  trom  Pierrehuirv 
beri  3l  Widnall  (I0|  while  ihe  dots  are  the  present  simulations. 
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Abstract 

The  response  of  circular  air  jets  to  combined  stream- 
wise  and  transverse  acoustic  excitation  is  described. 
The  jet  evolution  and  structure  are  documented  by 
flow  visualisation.  When  the  ratio  of  axial  to  trans¬ 
verse  excitation  frequencies  is  two,  the  jet  'bifurcates* 
into  two  distinct  jets.  The  spreading  angle  increases 
dramatically  with  excitation  amplitude. 

Nomenclature 

A  excitation  amplitude,  p' /ZqU 
D  noisle  exit  diameter 

/  excitation  frequency 

L  sound  pressure  level 

p'  acoustic  pressure 

R  fa/ft 

Re  Reynolds  number,  UD/v 
St  Strouhal  number, 

U  mean  exit  velocity 

Z  acoustic  impedance 

v  kinematic  viscosity 

Subscripts 
a  axial 

o  ambient 

i  transverse 

Introduction 

The  shear  layer  of  a  jet  naturally  rolls  up  into  dis¬ 
tinct  vortex  rings  due  to  the  Kelvin-Helmholts  insta¬ 
bility.  These  rings  combine  to  form  larger  structures 
which  also  pair  or  breakdown.  This  process  of  vortex 
formation  and  pairing  is  somewhat  random  in  the  un¬ 
excited  jet.  However,  by  proper  axial  excitation  one 
can  control  the  frequency  at  which  the  vortices  form 
and  subsequently  pair.  This  is  true  for  both  the  thin, 
closely-spaced  rings  near  the  exit  as  well  as  for  the 
larger  rings.1  Furthermore,  one  can  use  this  single¬ 
mode  forcing  to  control  jet  growth  by  enhancing  or 
suppressing  vortex  pairing.2 
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This  vortex  pairing  process  has  also  been  studied  in 
the  mixing  layer,  where  it  was  found  that  forcing  at 
a  subharmonic  of  the  most  amplified  frequency  results 
in  the  amalgamation  of  several  vortices.3  Only  very 
low  levels  of  forcing  are  required  to  cause  this  interac¬ 
tion.  This  pairing  mechanism  is  known  as  a  'collective 
interaction* . 

The  study  of  excited  jets  has  not  been  limited  to 
single-  frequency  plane  wave  excitation.  A  study  of 
active  cancellation  of  pure  tones  in  jets  involved  simul¬ 
taneous  perturbations  at  two  frequencies.4  The  ampli¬ 
fication  of  higher  order  modes  and  their  contribution 
to  broadband  jet  noise  has  also  been  considered.5  8  7 
Other  studies  of  excited  jets  have  focused  on  high¬ 
speed  jets*  and  jets  issuing  from  asymmetric  noxtles®. 

Properly  combining  axial  and  asimuthal  modes  can 
dramatically  alter  the  structure  and  momentum  trans¬ 
port  of  round  jets.10  This  type  of  dual-mode  forcing 
can  cause  a  jet  to  split  into  two  distinct  streams  (bi¬ 
furcating  jet)  or  to  explode  into  a  shower  of  vortex 
rings  (blooming  jet).  Bifurcating  jets  occur  when  the 
ratio  (R)  of  the  axial  to  transverse  frequency  is  two. 
When  R  is  a  non-integer  between  1.6  and  3.2,  blooming 
jets  are  observed.  These  phenomena  occur  only  within 
a  small  range  of  Strouhal  numbers,  and  within  that 
range  the  spreading  angle  increases  with  St.  These 
dramatic  changes  in  jet  development  were  discovered 
through  a  previous  experiment  in  water  at  a  Reynolds 
number  of  4,000. 11  A  vortex-filament  simulation  of  bi¬ 
furcating  jets  suggested  the  Strouhal  number  effect 
and  identified  vortex  interactions  as  a  key  mechanism 
in  this  flow.13 

Potential  applications  Athis  dual-mode  excitation 
involve  air  jets  at  higher  Reynolds  numbers.  Thus, 
this  work  extends  the  study  of  bifurcating  and  bloom¬ 
ing  jets  to  air  flows.  The  objective  here  is  to  control 
jet  flows  through  properly  organising  the  vortex  struc¬ 
tures  and  to  study  the  influence  of  excitation  ampli¬ 
tudes. 

This  paper  focuses  specifically  on  bifurcating  jets 
at  Re  of  10,000  and  20,000.  The  acoustic  excitation 
system  and  flow  visualisation  technique  are  described. 
Finally,  some  preliminary  results  on  the  influence  of 
excitation  amplitudes  are  discussed. 

Experimental  Apparatus  and  Approach 

Flow  System 

The  experimental  apparatus  is  shown  schematically 
in  Fig.  1.  Air  flows  through  a  porous  bronse  cylin- 
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der  into  the  plenum  and  is  directed  through  a  honey¬ 
comb  dink  upstream  of  the  nonle.  A  plexiglas  coupling 
rigidly  attaches  the  nosile  to  the  plenum.  The  noss’e 
exit  diameter  is  2.15  cm,  and  the  area  contraction  ra¬ 
tio  is  25  to  1.  The  mean  exit  velocities  are  7  and  14 
m/s,  which  correspond  to  Reynolds  numbers  of  10,000 
and  20,000,  respectively.  The  streamwise  turbulence 
intensity,  measured  with  a  hot-wire  on  the  centerline 
of  the  exit  plane,  is  1 %  of  the  mean  exit  velocity.  The 
shear  layer  thickness,  measured  with  a  hot-wire  0.1  di¬ 
ameter  from  the  exit,  is  3%  of  the  diameter  at  Re  of 
10,000. 

Flow  Visualization 

The  plexiglas  coupling  shown  in  Fig.  1  also  serves 
as  the  passage  for  smoke  injection.  Cigar  smoke  en¬ 
ters  through  four  ports  into  the  coupling  and  forms  a 
thin  cylindrical  sheet  as  it  exits  from  the  annular  slot 
in  a  direction  tangential  to  the  core  flow.  Alternately, 
smoke  can  be  channeled  through  only  two  opposite 
ports  to  provide  an  approximate  cross-section  of  the 
jet.  By  injecting  smoke  into  the  boundary  layer,  in¬ 
stead  of  seeding  the  entire  jet,  one  clearly  marks  the 
shear  layer  and  vortex  structures. 

Both  time-averaged  and  instantaneous  pictures  of 
the  flow  are  taken.  Time-averaged  pictures  are  ob¬ 
tained  by  strobe  illumination.  The  strobe  is  triggered 
at  a  submultiple  of  the  excitation  frequency  and  is  thus 
able  to  phase-lock  the  evolution  of  the  vortex  struc¬ 
tures.  The  images  are  recorded  by  a  video  camera.  A 
10-watt,  copper-vapor,  pulsed  laser  is  focused  into  a 
thin  sheet  to  provide  instantaneous  cross-sections  of 
the  jet.  The  laser  is  triggered  to  provide  one  pulse  for 
each  exposure  taken  with  a  35mm  camera  and  stan¬ 
dard  print  film  (ASA  400).  Since  the  pulse  duration  is 
only  30  ns,  this  technique  provides  detailed  images  of 
even  high  speed  flows. 

Acoustic  Excitation 

The  bifurcating  jet  experiment  in  water  used  a  me¬ 
chanical  excitation  system  operating  at  frequencies 
around  lOHs.11  A  diaphragm  on  the  piston  driving  the 
flow  produced  the  axial  perturbation.  The  asimuthal 
mode  was  introduced  by  moving  the  nossle  tip  in  a 
small  circular  orbit  about  the  nominal  jet  centerline. 
From  the  definitions  of  the  Reynolds  and  Strouhal 
numbers,  the  excitation  frequency  can  be  written  as 
/<,  =  vReSt/D2.  This  shows  that  switching  from  wa¬ 
ter  to  air,  while  keeping  Re,  St,  and  D  constant,  re¬ 
quires  an  order  of  magnitude  increase  in  /„. 

Since  the  frequencies  required  in  this  experiment  are 
several  hundred  herts,  an  acoustic  excitation  system 
was  developed.  Loudspeakers  produce  both  the  ax¬ 
ial  and  transverse  excitations.  A  speaker  mounted 
on  the  bottom  of  the  plenum  generates  the  stream- 
wise  excitation.  Four  external  speakers  (Morel,  15- 
cm  woofers,  150W)  surrounding  the  nossle  provide  the 
cross-stream  excitation.  By  mounting  these  speakers 


in  the  plane  of  the  nossle  exit,  one  presents  a  flat  sur¬ 
face  to  the  entrained  flow. 

The  input  signals  of  the  external  speakers  are  sine 
waves  of  equal  frequency  and  amplitude  but  different 
phase.  The  signal  of  a  speaker  differs  from  that  of 
an  adjacent  one  by  90°.  The  signal  of  the  opposite 
speaker  is  180°  out  of  phase.  The  input  of  the  internal 
speaker  is  a  sine  wave  of  a  different  frequency,  and  its 
phase  relative  to  the  set  of  four  speakers  is  variable. 

One  pair  of  opposite  external  speakers  is  sufficient 
to  generate  the  transverse  forcing  required  to  produce 
the  bifurcating  jet.  However,  the  jet  will  split  only  if 
the  phase  shift  between  the  pair  of  external  speakers 
and  the  internal  one  is  properly  set.  The  peaks  of 
the  transverse  excitation  must  coincide  with  the  peaks 
of  the  axial  forcing.  With  all  four  external  speakers 
on,  the  phase  shift  simply  deteimines  the  diametrical 
plane  in  which  the  jet  bifurcates.  The  results  reported 
here  are  produced  with  only  the  internal  speaker  and 
one  external  pair  of  opposite  speakers  properly  phased. 
No  signal  is  fed  into  the  other  pair. 

The  sound  pressure  levels  are  measured  with  a  2.4 
cm  microphone  centered  one  diameter  from  the  jet  exit 
and  oriented  perpendicular  to  the  plane  containing  the 
speaker  pair.  The  purpose  of  this  measurement  is  to 
specify  the  excitation  amplitudes  and  not  to  obtain 
the  acoustic  response  of  the  excited  jet.  In  order  to 
determine  the  forcing  amplitude  prior  to  amplification 
by  jet  instability  mechanisms,  all  acoustic  measure¬ 
ments  are  made  with  no  flow.  The  two  modes  of  exci¬ 
tation  are  measured  separately.  The  transverse  excita¬ 
tion  amplitudes  that  are  reported  are  those  resulting 
from  both  speakers.  The  corresponding  transverse  am¬ 
plitude  from  a  single  speaker  is  typically  15  dB  higher 
than  that  of  the  pair  due  to  the  absence  of  cancellation 
by  the  other  speaker. 

The  spreading  angles  are  determined  from  the  flow 
visualisation.  When  instantaneous  pictures  are  used, 
several  realisations  are  compared  to  determine  the  an¬ 
gle.  The  spreading  angle  of  the  bifurcating  jet  is  de¬ 
fined  as  the  angle  formed  by  the  outer  edges  of  the  two 
branches  of  the  jet. 

Renulta  and  Discussion 

The  results  discussed  here  are  of  a  preliminary  na¬ 
ture,  yet  they  provide  valuable  insight  into  the  struc¬ 
ture  and  amplitude  dependence  of  bifurcating  jets. 
The  Strouhal  number,  defined  as  faD/U,  is  around  0.6 
in  all  these  cases.  This  value  of  St  was  found  in  pre¬ 
vious  work  to  produce  the  largest  spreading  angles.11 
The  non-dimensional  amplitudes,  Aa  and  At,  are  de¬ 
fined  as  the  fluctuating  acoustic  velocity  divided  by 
the  mean  jet  velocity.  The  acoustic  velocity  is  esti¬ 
mated  by  dividing  the  measured  acoustic  pressure  by 
the  characteristic  impedance  of  air.  Thus,  Aa  and  At 
are  defined  as  p’a/ZoU  and  p't/ZoU,  respectively. 


With  both  excitation*  off,  we  simply  observe  the 
natural  jet  whose  shear  layer  rolls  up  into  vortices  in  an 
irregular  manner  (Fig.  2).  Turning  the  axial  excitation 
on  locks  the  vortex  roll-up  to  the  forcing  frequency  as 
seen  in  Fig.  3. 

Cross-sections  through  the  center  of  the  jet  show 
interesting  differences  in  jet  response  to  Aa.  At  low 
levels  of  streamwise  forcing,  one  controls  the  vortex 
formation  frequency  by  causing  a  periodic  collective 
interaction  of  the  thin  closely-  spaced  rings  near  the  jet 
exit.  The  axially-pulsed  jet  at  Aa  —  0.03%  (La  =  99 
dB)  is  an  example  of  this  (Fig.  4).  At  higher  levels, 
Aa  =  0.12%  ( La  =  105  dB),  the  thin  vortex  rings  are 
no  longer  visible,  and  instead,  one  sees  the  initial  for¬ 
mation  of  rings  whose  spacing  is  on  the  order  of  the 
jet  diameter  (Fig.  5).  The  stronger  excitation  causes 
a  tight  roll-up  of  the  shear  layer.  The  structure  of 
the  vortex  core  shows  many  interfaces  between  the  jet 
fluid  and  the  entrained  fluid.  Several  diameters  from 
the  exit,  vortex  cores  exhibit  deformation  and  elonga¬ 
tion  due  to  the  pairing  process.  The  same  amplitude 
dependence  is  seen  in  jets  having  the  same  Reynolds 
number,  in  which  case  the  thickness  of  the  shear  layer 
is  not  an  issue. 

The  bifurcating  jet  has  a  strikingly  different  struc¬ 
ture.  Figure  6  displays  a  cross-section  of  a  bifurcat¬ 
ing  jet  in  the  plane  of  bifurcation  at  Re  =  10,000 
and  St  =  0.55.  The  jet  appears  similar  to  an  axially- 
excited  jet  near  the  exit  except  for  the  slight  displace¬ 
ment  and  tilt  of  the  vortices.  Farther  downstream, 
however,  the  flow  abruptly  splits  into  two  distinct  jets 
with  an  included  angle  of  70°.  Initially  adjacent  vor¬ 
tex  rings  propagate  along  different  branches  of  the  jet. 
This  results  in  the  jet  fluid  being  stretched  back  and 
forth  between  the  two  branches  of  the  jet.  Similar  be¬ 
havior  is  seen  at  Re  =  20,000  (Fig.  7).  Since  the 
Reynolds  number  is  higher  and  amplitudes  are  lower, 
the  vortex  cores  are  not  aa  distinctly  organised  and 
the  jet  becomes  fully  turbulent  sooner. 

Changes  in  spreading  angle  due  to  different  trans¬ 
verse  amplitudes  are  seen  in  Fig.  8.  At  low  amplitudes 
the  vortices  appear  to  be  displaced  laterally  rather 
than  distributed  on  two  separate  branches.  The  angles 
corresponding  to  different  levels  of  axial  and  transverse 
forcing  are  compared  in  Fig.  9.  Increasing  At  while 
keeping  Aa  fixed  causes  the  spreading  angle  to  increase 
from  24°  up  to  52°.  Increasing  Aa  results  in  angles  as 
high  as  70°. 

Though  more  data  would  be  required  to  form  an 
empirical  relation  between  angle  and  amplitude,  these 
results  provide  the  first  insights  into  the  nature  of  that 
dependence.  The  spreading  angle  clearly  increases 
with  both  amplitudes.  One  would  intuitively  expect 
the  angle  to  increase  with  At  since  that  perturbation 
is  aligned  with  the  direction  of  bifurcation.  The  fact 
that  the  angle  dramatically  increases  with  A„  would 
be  harder  to  anticipate.  It  appears  that  higher  values 


of  Aa  concentrate  the  vorticity  in  the  shear  layer  more 
strongly.  Thus,  the  vortices  are  less  diffuse  and  prob¬ 
ably  have  a  higher  circulation.  This,  in  turn,  amplifies 
the  vortex  interactions,  resulting  in  wider  spreading 
angles.  One  would  expect  to  reach  a  saturation  level 
beyond  which  the  angle  no  longer  increases  with  either 
amplitude.  However,  that  level  did  not  appear  to  be 
reached  in  these  experiments. 

These  results  confirm  and  expand  the  previous  un¬ 
derstanding  of  the  governing  mechanisms  in  this  flow. 
The  bifurcation  phenomenon  has  been  regarded  as  the 
result  of  vortex  interactions.11,12  Eccentric  rings  mu¬ 
tually  induce  each  other  to  tilt  away  from  the  center- 
line.  This  process  amplifies  until  the  jet  becomes  two 
separate  trains  of  collinear  rings.  The  present  results 
are  in  agreement  with  that  model.  Additionally,  this 
experiment  demonstrates  that  the  vortex  interactions 
can  be  intensified  by  increasing  either  excitation  am¬ 
plitude. 

Conclusions 

An  acoustic  excitation  system  capable  of  generating 
both  streamwise  and  cross-stream  perturbations  was 
developed  for  round  air  jets.  The  effect  of  this  dual¬ 
mode  forcing  at  different  Reynolds  numbers  and  ex¬ 
citation  amplitudes  is  visually  documented  These  re¬ 
sults  demonstrate  that  dual-mode  excitation  can  cause 
a  jet  to  split  into  two  distinct  jets  with  a  spreading  an¬ 
gle  as  high  as  70°.  The  Reynolds  number  range  for  bi¬ 
furcation  has  been  extended  to  as  high  as  20,000  with 
excitation  frequencies  up  to  386  Hs. 

This  experiment  agrees  with  previous  understand¬ 
ing  of  the  bifurcating  jet  structure  and  past  models  of 
the  vortex  interactions.  Additionally,  this  work  shows 
that  the  spreading  angle  of  the  bifurcating  jet  increases 
with  both  axial  and  transverse  amplitudes. 
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(a) 


Transverse  Excitation  Amplitude 


Figure  9.  Variation  of  spreading  angle  with  excitation 
amplitude. 
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(b) 


Figure  8.  Video  images  of  bifurcating  jets  at  different 
transverse  excitation  amplitudes. 

Re  =  10,000,  St  =  0.60,  fa  =  192  Hz, 
ft  =  96  Hz,  and  Aa  =  0.05%. 

(a)  At  =  0.01%,  (b)  At  =  0.02%,  (c)  At  =  0.04%. 
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ABSTRACT 

Dual-mode,  dual-frequency  acoustic  excitation  of 
round  air  jets  is  described.  The  jet  evolution  and 
structure  is  documented  by  flow  visualization  at 
velocities  up  to  75  m/s  and  Reynolds  numbers  to 
100,000.  The  ratio  of  the  axial  to  helical  excitation 
frequencies  is  exactly  two.  This  type  of  forcing  causes 
the  jet  to  spread  dramatically  in  one  plane.  The 
spreading  angle  increases  with  excitation  amplitude  to 
angles  as  high  as  70  degrees. 


INTRODUCTION 

The  sensitivity  of  jets  to  sound  has  fascinated 
researchers  for  many  decades.  Brown  demonstrated 
that  laminar  jets  develop  vortex  structures  and  increase 
in  spreading  angle  in  response  to  acoustic  excitation  at 
various  critical  frequencies  (1).  Most  current  work 
focuses  on  jets  consisting  of  an  axisymmetric  shear 
layer  and  a  potential  core.  In  this  case  the  shear  layer 
rolls  up  to  form  distinct  vortex  rings  unlike  the  laminar 
jets  of  Brown  which  form  vortex  structures  comprised 
of  the  entire  jet  stream. 

The  vortex-formation  frequency  can  be  fixed  by 
axial  excitation.  This  single-mode  forcing  can  control 
jet  growth  by  enhancing  or  suppressing  vortex  pairing 
(2,3).  The  "collective  interaction"  of  several  vortices 
due  to  subharmonic  forcing  has  been  demonstrated  in 
mixing  layers  (4).  This  type  of  forcing  has  also  been 
applied  to  high  speed  jets  (5)  and  jets  issuing  from 
asymmetric  nozzles  (6). 

The  study  of  excited  jets  has  involved  multiple- 
frequencies  and  multiple-modes.  A  study  of  active 
cancellation  of  pure  tones  in  jets  involved  simultaneous 
perturbations  at  two  frequencies  (7).  Other  work  has 
focused  on  the  changes  in  initial  shear  layer 
development  in  response  to  multiple-frequency  forcing 
(8). 

Properly  combining  axial  and  first-order  helical 
modes  can  dramatically  alter  the  structure  and 
momentum  transport  of  round  jets.  This  type  of  dual- 
mode  forcing  can  cause  a  jet  to  split  into  two  distinct 
streams  when  the  ratio  of  the  axial  to  helical  frequency 
is  two.  These  phenomena  occur  only  within  a  small 


range  of  Strouhal  numbers,  and  within  that  range  the 
spreading  angle  increases  with  St.  These  dramauc 
changes  in  jet  development  were  discovered  in  a  study 
of  mechanically  perturbed  round  water  jets  at  a 
Reynolds  number  of  4,000  (9).  Similar  jet 
development  was  observed  in  acoustically-excited  air 
jets  in  the  Reynolds  number  range  of  10,000  to  20,000 
(10). 

This  work  extends  the  study  of  bifurcating  air  jets  to 
Re  of  100,000.  Since  the  apparatus  described  in  Ref. 
10  is  inadequate  to  produce  the  high  levels  of 
excitation  required  at  higher  velocities,  a  new  acoustic 
excitation  system  was  developed.  The  response  of  the 
jet  at  various  Reynolds  numbers  is  visually 
documented. 


EXPERIMENTAL  APPARATUS  AND  APPROACH 


Flow  System 

The  experimental  apparatus  is  shown  schematically 
in  Fig.  I.  Air  flows  through  a  porous  bronze  cylinder 
into  the  plenum  and  out  through  a  2-cm-diametcr 
nozzle.  The  two-piece  nozzle  has  a  carefully  machined 
fifth-order  polynomial  profile  with  zero  slope  at  inlet 
and  exit.  The  transition  from  the  lower  to  the  upper 
half  of  the  nozzle  occurs  at  the  inflection  point  of  the 
profile.  The  area  contraction  ratio  is  25  to  1 ,  and  the 
length-to-diameter  ratio  is  5. 

The  jet  exit  is  positioned  in  the  center  of  a  60  cm  x 
60  cm  fiat  panel  and  is  flush  with  this  panel.  Baffles 
made  of  particle  board  and  Sonex  acoustical  foam 
surround  the  jet  on  all  four  sides.  These  baffles  are 
located  1  m  from  the  jet  on  each  side.  The  fume  hood 
located  1.3  m  above  the  jet  is  lined  with  Sonex  foam. 
Entrained  air  flows  into  the  test  cell  through  the  0.8-m 
gap  between  the  baffles  and  the  floor. 

Flow  Visualization 

A  small  0.5-mm  annular  gap  between  the  two 
halves  of  the  nozzle  provides  a  passage  for  injecting  a 
fluid  marker  into  the  boundary  layer.  In  these 


in  smoke  reaching  the  jet  centerline  about  one  diameter 
closer  to  the  jet  exit  (Fig.  8). 

An  axially-excited  jet  at  Re  of  100,000  is  shown  in 
Fig.  9.  The  high  level  of  forcing  causes  the  shear  layer 
to  turn  turbulent  at  the  exit.  At  low  levels  of  excitation 
(below  120  dB),  the  shear  layer  remains  laminar  for  the 
same  distance  as  in  the  unforced  jet  (Fig.  6).  As 
previously  discovered  (10).  the  spreading  angle 
increases  with  excitation  amplitude  to  angles  as  high  as 
70  degrees  (Figs.  10  and  11).  At  the  lower  Reynolds 
numbers,  increasing  the  helical  excitation  amplitude 
beyond  a  certain  level  does  not  increase  the  spreading 
angle  and  in  some  instances  seems  to  reduce  it.  This  is 
true  regardless  of  how  the  axial  excitation  is 
introduced.  The  excitation  levels  required  to  reach  this 
saturation  level  is  found  to  increase  with  Reynolds 
numbers. 

Adjusting  the  phase  between  the  axial  and  helical 
signals  rotates  the  plane  in  which  the  jet  bifurcates.  By 
rotating  the  jet  such  that  it  bifurcates  in  the  plane 
perpendicular  to  the  light  sheet,  one  obtains  the  image 
in  Fig.  12.  Since  the  jet  fluid  moves  away  from  this 
bisecting  plane,  the  smoke  density  along  the  jet 
centerline  decreases  rapidly.  The  striking  difference 
between  the  cross-sections  in  Figs.  10  and  1 1 
demonstrates  the  fact  that  the  bifurcating  jet  spreads 
rapidly  in  one  plane  rather  than  axisymmetrically. 


CONCLUSIONS 

The  effect  of  combined  axial  and  helical  excitations 
on  the  development  of  round,  turbulent  air  jets  has 
been  studied  by  flow  visualization.  Setting  the  axial 
frequency  at  exactly  double  the  helical  frequency 
causes  the  jet  to  spread  rapidly  in  one  diametrical 
plane.  Increasing  the  amplitude  of  the  helical 
excitation  results  in  spreading  angles  as  high  as  70 
degrees.  However,  increases  above  a  certain  level  do 
not  further  amplify  the  spreading  angle.  Higher  levels 
of  excitation  are  required  to  reach  this  saturation  level 
at  higher  Reynolds  numbers.  This  work  demonstrates 
that  dual-mode  forcing  can  cause  jets  to  bifurcate  at 
velocities  up  to  75  m/s  and  Reynolds  numbers  up  to 
100,000. 
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1.  Page  3:  "Figs.  10  and  1 1  demonstrates '  should 
be  "Figs.  10  and  12  demonstrates. 


2.  Page  4:  "Figure  4 - (N  =  32)"  should  be 

"Figure  4.  .  .  .  (N  =  8).'" 
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